Is there Correlation between Fine Structure and Dark Energy Cosmic Dipoles? 
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We present a detailed analysis (including redshift tomography) of the cosmic dipoles in the 
Keck+VLT quasar absorber and in the Union2 Snla samples. We show that the fine structure 
constant cosmic dipole obtained through the Keck+VLT quasar absorber sample at 4.1a level is 
anomalously aligned with the corresponding dark energy dipole obtained through the Union2 sam- 
ple at 2a level. The angular separation between the two dipole directions is 11.3° ± 11.8°. We 
use Monte Carlo simulations to find the probability of obtaining the observed dipole magnitudes 
with the observed alignment, in the context of an isotropic cosmological model with no correla- 
tion between dark energy and fine structure constant a. We find that this probability is less than 
one part in 10 6 . We propose a simple physical model (extended topological quintessence) which 
naturally predicts a spherical inhomogeneous distribution for both dark energy density and fine 
structure constant values. The model is based on the existence of a recently formed giant global 
monopole with Hubble scale core which also couples non-minimally to electromagnetism. Aligned 
dipole anisotropies would naturally emerge for an off-centre observer for both the fine structure 
constant and for dark energy density. This model smoothly reduces to ACDM for proper limits of 
its parameters. Two predictions of this model are (a) a correlation between the existence of strong 
cosmic electromagnetic fields and the value of a and (b) the existence of a dark flow on Hubble 
scales due to the repulsive gravity of the global defect core ('Great Repulser') aligned with the dark 
energy and a dipoles. The direction of the dark flow is predicted to be towards the spatial region of 
lower accelerating expansion. Existing data about the dark flow are consistent with this prediction. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



> 
in 
in 
o 
^t- 

\6 
o 



X 



1. INTRODUCTION 

According to the cosmological principle, the Universe 
is homogeneous and isotropic on scales larger than a few 
hundred Mpc. The main source of evidence which sup- 
ports this assumption comes from the Cosmic Microwave 
Background (CMB) which appears to be isotropic to a 
high degree up to a dipole term which is assumed to be 
due to our motion with respect to the CMB frame. How- 
ever, there has been some recent observational evidence 
which could be interpreted as a hint for deviations from 
large scale statistical isotropy. Such evidence includes 
alignment of low multipoles in the CMB angular power 
spectrum [l|, large scale velocity flows @, y| and large 
scale alignment in the QSO optical polarisation data [4[ 
(see Ref. Q for an interesting related theoretical model) . 
These effects appear to persist on scales of 1 Gpc or larger 
and could constitute early hints for a deviation from the 
FLRW metric on large cosmological scales and the exis- 
tence of a cosmological preferred axis. This possibility is 
further enhanced by the fact that the anisotropy direc- 
tions implied by these observations appear to be abnor- 
mally close to each other Q . 
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The above hints for cosmological anisotropy have moti- 
vated searches for deviations from the cosmological prin- 
ciple by considering the angular distribution of luminos- 
ity distances of Type la supernovae (Snla) in the red- 
shift range z € [0.015, 1.4] [l-[l4]|. Even though all these 
studies are consistent with isotropy, in many of them, 
a mild evidence_(l(7 — 2a) of anisotropic expansion was 
found [H, \l4 1 mainly coming from low redshift data, 
while in others [9 , [111, [n| no evidence of anisotropy was 
found. The inability of the later studies to pick up any 
anisotropy is perhaps due to the methods and data used 
which were not sensitive enough to particular types of 
anisotropy. 

Additional hints for such possible deviations from the 
cosmological principle have recently been obtained by the 
angular distribution of the fine structure constant a in 
the redshift range z £ [0.2223, 4.1798] as measured by the 
quasar absorption line spectra using the many multiplet 
method [l5j. If in the case of Snla the dipole anisotropy 
was mild (about 1 — 2a), in the case of the fine structure 
constant the anisotropy has been found to be significant 
(4.1a). 

Some earlier studies had claimed possible variation of 
the fine structure 'constant' with time Hi . This possibil- 
ity has led to extensive theoretical modelling in the litera- 
ture so far [It], [3 with emphasis on the possible connec- 
tion of this variation with dark energy (quintessence) fl9| . 
However, there has been comparatively less interest in 
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the possibility of spatial variation of a (see however 
[20I l2lj for recent studies) and its connection with dark 
energy. 

The anisotropy analysis of Ref. [6| fo r the Snla sample 
was based on the Union2 dataset [23| which consists of 
557 Snla. A hemisphere comparison method was used to 
find the hemisphere pair with maximal anisotropy with 
respect to ACDM fits. The maximum anisotropy direc- 
tion was found to be towards (/, b) = (309°, 18°) but the 
magnitude of this dark energy anisotropy was found to 
be consistent with statistical isotropy at the 2a level. 
Similar results were obtained in Ref. [7j where a redshift 
tomography also revealed that most of the contribution 
to the mild dark energy dipole comes from the low red- 
shift Snla. 

The anisotropy analysis of the fine structure constant 
a [IH is based on a large sample of quasar absorption-line 
spectra (295 spectra) obtained using UVES (the Ultravi- 
olet and Visual Echelle Spectrograph) on the VLT (Very 
Large Telescope) in Chile and also previous observations 
at the Keck Observatory in Hawaii. An apparent varia- 
tion of a across the sky was found. It was shown to be 
well fit by an angular dipole model (— ) = A cos 6 + B 
where 9 is the angle with respect to a preferred axis and 
A, B are the dipole magnitude and an isotropic monopole 
term. The dipole axis was found to point in the direc- 
tion (/, b) = (331°, —14°) and the dipole amplitude A was 
found to be A = (0.97 ±0.21) x 10" 5 . The statistical sig- 
nificance over an isotropic model was found to be at the 
4.1(7 level. The analysis of Ref. [15| has received criti- 
cism [22[ based mainly on the fact that its quasar sample 
combines two datasets (Keck and VLT) with different 
systematic errors which have a small overlapping subset 
and cover opposite hemispheres on the sky. The axis con- 
necting these two hemispheres has similar direction with 
the direction of the obtained dipole. The response of the 
authors of Ref. was based on the fact that in the 
equatorial region of the dipole, where both the Keck and 
VLT samples contribute a number of absorbers, there is 
no evidence for inconsistency between Keck and VLT. 

The controversy about the possible problems in the 
analysis of Ref. [15| and the angular proximity between 
the dark energy axis of Ref. [6[ and the (^) axis of 
Ref. 15[ constitutes the motivation to analyse both the 
Snla and the quasar datasets in a similar and consis- 
tent manner. Thus we re-analysed both datasets and fit 
them to the same dipole+monopole ansatz of the form 
A cos 8 + B. This type of anisotropy fit is different from 
the corresponding Snla fits of previous studies. Our goal 
is to address the following questions: 

1. What are the best fit dipoles (magnitudes A and 
directions in galactic coordinates) for the Union2 
and Keck+VLT samples? What is the angle be- 
tween the two dipole directions? 

2. How likely is it to obtain these dipole magnitudes in 
the context of an isotropic underlying model? How 
likely is it to obtain the observed angle between the 



dipoles if the two underlying models were isotropic 
and uncorrelated? We address these questions by 
producing a large number of Monte Carlo isotropic 
datasets simulating the Union2 and the Keck+VLT 
samples under the assumption of isotropic and un- 
correlated underlying models. We then compare 
the obtained probability distributions for the dipole 
magnitudes and angles with the observed magni- 
tudes and angle. 

3. How do the answers to the above questions change 
if we consider three different redshift slices (bins) 
for each dataset (low, medium and high redshift) 
with approximately equal number of datapoints in 
each bin? Is there a particular redshift range where 
the dark energy and the fine structure dipoles are 
more prominent and how is this range related with 
the quality of the data in each bin? 

These questions are addressed in detail in the follow- 
ing sections. In particular, the structure of this paper 
is the following: in the next section we derive the mag- 
nitudes and directions of the best fit dark energy and 
fine structure dipoles for the full Union2 and Keck+VLT 
datasets thus addressing the above question 1. We also 
perform 10 4 Monte Carlo simulations of the Union2 and 
Keck+VLT datasets based on an isotropic best fit ACDM 
model and on an isotropic best fit monopole model re- 
spectively. We then use these simulations to address the 
above question 2. In section [3] we perform a redshift 
tomography to address question 3 and find the redshift 
range where the dipoles appear to be more pronounced. 
In section|3]we discuss a physical model that could repro- 
duce the observed dipole alignment. Finally, in section [5] 
we conclude, summarise our basic results and discuss fu- 
ture prospects of the present work. 



2. COSMIC DIPOLES: DATA ANALYSIS AND 
MONTE CARLO SIMULATIONS 

A. Fine Structure Constant Dipole 

The full Keck+VLT sample consists of 295 quasar 
absorption line spectra in a redshift range z € 
[0.2223, 4.1798]. It has been analysed in detail in Ref. [3 
where the redshift of each absorber is presented along 
with fine structure constant deviation (—) = az ~ a ° 
where a z is the value of a measured at redshift z using 
the many multiplet method (lq | and «o is the value of a 
measured in the laboratory. The positions of the quasars 
in equatorial coordinates are also presented. 

In order to fit the Keck+VLT dataset to a dipole 
anisotropy we proceed as follows: 

• We convert the equatorial coordinates of each 
quasar to galactic coordinates. 

• We find the Cartesian coordinates of the unit vec- 
tors hi corresponding to each quasar with galactic 
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coordinates (l,b). We thus have 
hi = cos(6i) cos(li)i + cos(6i) sin(Zj)j + sin(&i)fc (2.1) 

• We use the dipole+monopole angular distribution 
model 
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A cos 6 + B 



(2.2) 



where cos 8 is the angle with the dipole axis defined 
by the vector 



such that 



D = c\i + c 2 j + c 3 k 



■ D = A cos t 



(2.3) 



(2.4) 



We fit the Keck+VLT dataset to a dipole 
anisotropy model (|2.2j) using the maximum likeli- 
hood method i.e. minimising 



X 2 (D,B) 
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Acose.-B]' 



(2.5) 



where (—). and cr, are obtained from the 
r - hi 

Keck+VLT dataset [15] and ay an <2 is an internal 
random error, assumed to be the same for all data 
points and representing an estimate of the aggre- 
gation of all additional random errors. We fix the 
value of (Jrcmd by requiring that at the best fit 
X 2 (D, B) per degree of freedom is about unity. The 
required value of a ran d is 1.0 x 10~ 5 in agreement 
with the corresponding value used in Ref. fl5| . 

• The magnitude and direction of the best fit dipole 
in galactic coordinates is obtained from the best 
fit a coordinates (e.g. A = c\ + c\ + cf) and the 
corresponding la errors are obtained using the co- 
variance matrix approach. 

Our result for the best fit dipole direction and mag- 
nitude is consistent with the corresponding results of 
Ref. QJ] (see also Table EJ). We find A fs = (1.02±0.25) x 
10" 5 with direction (b = -11.7°±7.5°, I = 320.5°±11.8°) 
while for the best fit monopole term we have Bf s = 
(— 2.2± 1.0) x 10~ 6 . This result shows that the isotropic 
model A = is more than 4cr away from the best fit 
value of the dipole magnitude. The Keck+VLT dataset 
along with the best fit direction of the dipole in galactic 
coordinates is shown in Fig. [I] By the definition of (— ) 
and the construction of the dipole model, the obtained 
dipole direction shown in Fig. [1] is the direction towards 
larger values of the fine structure constant a. 

In an effort to better analyse the above results for the 
best fit dipole and its errors we have constructed a Monte 
Carlo simulation obtained from the Keck+VLT dataset 
under the assumption of an isotropic monopole model. 



Such a simulation is aimed at providing the probability 
distribution of the dipole magnitude and direction un- 
der the assumption of an isotropic monopole model and 
through that, the probability of obtaining the actually 
measured values in the context of an isotropic model. In 
order to construct the Monte Carlo simulation we pro- 
ceed as follows: 

• We define a Gaussian random selection function 
g(n,a) which returns a random number from a 
Gaussian probability distribution with mean /j, and 
variance a 1 . 

• We fit the Keck+VLT dataset to an isotropic 
monopole model obtained from eq. (|2.2j) by setting 
A = 0. We find for the best fit monopole term: 
Bf s - m = (-0.19 + 0.10) x 10~ 5 . 

• We construct the isotropic Monte Carlo version of 
the Keck+VLT dataset by keeping fixed the di- 
rection of each quasar and assigning to each ab- 
sorber an isotropic randomised fine structure con- 
stant variation obtained as: 



Aa 



AIC 



g(Bf s ^ mi <Ti) + g(0,a ran d) (2.6) 



• We construct 10 4 such Monte Carlo datasets and 
obtain the probability distribution of the dipole 
magnitude as well as the corresponding dipole di- 
rections. We thus find the number of isotropic 
datasets that have a dipole magnitude larger than 
the observed value of the dipole magnitude. 

The obtained probability distribution of the dipole mag- 
nitudes is shown in Fig. [5] along with the actually ob- 
served value of A. None of the 10 4 isotropic Monte Carlo 
datasets had a dipole magnitude as large as the one ob- 
served (or larger). We thus conclude that the proba- 
bility to obtain the observed dipole magnitude of the 
Keck+VLT dataset in the context of an isotropic model 
is less than 0.01% (3.9ct) in agreement with the covari- 
ance matrix error and with the result of Ref. [l5[ where 
the value 4.1cr was obtained. 



B. Dark Energy Dipole 

We perform a similar dipole+monopole fit using the 
Union2 data. Instead of (— ), which corresponds to fine 
structure constant deviations from its earth measured 
value, we use the distance modulus deviation from its 
best fit ACDM value 



Afj,(z)\ _ n{z) - fi(z) 



(2.7) 



where p, is the best fit distance modulus in the context of 
ACDM . The 557 Snla data points of the Union2 dataset 
are given in terms of the distance moduli 



fiobs(zi) = m obs (zi) - M 



(2. 
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(b=90) 
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FIG. 1: Keck+VLT datapoints and a-dipole direction. Datapoints in three different redshift bins are represented with different 
shapes. For comparison the direction of the Dark Energy dipole obtained from the best fit of the Union2 data is shown with a 
star. The light blue blob represents the 1-a error on the a-dipole direction. 



1000 
800 
600 
400 
200 



0.2 



0.4 



0.6 



0.8 



1.0 



dK/v00~ b ) 



FIG. 2: Distribution of a-dipole magnitudes obtained from 
the Monte Carlo simulation. The arrow points to the position 
of the observed best fit value for the a-dipole magnitude. 



where m b s is the apparent magnitude of each Snla and 
M is the absolute magnitude assumed to be common for 
all Snla after proper calibration. Assuming a ACDM 
parametrisation of the expansion rate 

H{zf = H$[n 0m (l + zf + (1 - Slom)] (2.9) 

the best fit distance modulus p>(z) is determined by min- 
imising 



2 co ^ \- [^QbsjZr) - Hth{Zj)} 2 
X l"0m,^0j — 2 



(2.10) 



i=i 



where a 1 ^ i are the distance modulus uncertainties which 
include both the observational and the intrinsic random 
magnitude scatter. The theoretical distance modulus is 
defined as 

l^th(zi) = m t h(zi) - M = 5log w (D L (z)) + ^ (2.11) 



where fj,o is a constant related to the Hubble parameter 
H = 100 h km/(sec • Mpc) by 



and 



fio = 42.38 — blogioh 



Dl(z) = (1 + z) f Z dz' H ° 

Jo H{z';ilom) 



(2.12) 



(2.13) 



is the Hubble free luminosity distance. A minimisation of 
X 2 (^om, A*o) using the Union2 dataset leads to the best fit 
parameter values flora = 0.269 ± 0.020 and fiQ = 43.16 ± 
0.01 which completely specify fi(Zi) and therefore 



An(zj) 



obs 



fi(zj) - n(zj) 



(2.14) 



for all Union2 datapoints. 

We now perform the same analysis as for the 
Keck+VLT data, replacing the (^) datapoints by the 

( A p(z) ) datapoints. In the Snla we set a r and = 
since the random intrinsic magnitude scatter has al- 
ready been included in the distance moduli errors <7j. 
We find the direction of the dark energy dipole to be 
(b = -15.1° ± 11.5°, I = 309.4° ± 18.0°). The magni- 
tudes of the dipole and monopolc terms are found to be 

A de = (1.3 ±0.6) x 10~ 3 (2.15) 
B de = (2.0 ±2.2) x 10~ 4 (2.16) 

The statistical significance of the dark energy dipole is at 
the 2cr level (significantly smaller that the 4cr of the fine 
structure constant dipole) but its direction is only 11° 
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away from the corresponding direction of the fine struc- 
ture constant dipole. The direction of the dark energy 

dipole along with the Union2 data ( A ^ z ^ ) are shown 

\ W Z »J J b s 

in Fig. |3] in galactic coordinates. The proximity of the 
two dipole directions is also made apparent in the same 
plot as well as by comparing with Fig. [TJ 

The direction of the dipole in Fig. [3] points towards 
brighter Snla compared to best fit isotropic ACDM . This 
implies less accelerating expansion in that direction only 
if H (which is related to /j,q in eq. (|2.1ip ) is assumed to 
be isotropic. This assumption was not made in Ref. @ 
where [1q was simultaneously fit along with fiom in each 
hemisphere. However, in the discussion of section 0] we 
will assume isotropic fiQ and therefore lower acceleration 
in the direction of brighter Snla. 

In an effort to determine the likelihood of the ob- 
served dark energy dipole magnitude combined with its 
angular proximity to the fine structure dipole we have 
performed a Monte Carlo simulation consisting of 10 4 
Union2 datasets constructed under the assumption of 
isotropic ACDM . Thus the distance modulus of point 
i is given by 

fJ,Mc(zi) = g(P>(zi),<ri) (2-17) 

where g is the Gaussian random selection function de- 
fined in the previous subsection and jl(zi) is the best 
fit distance modulus of the Union2 data in the context 
of ACDM at redshift Zi. It is thus straightforward to 

construct ( ) for each Monte Carlo dataset and 

V « z *) J MC 

obtain its best fit dipole direction and magnitude. In 
Fig.|4]we show the probability distribution of the dark en- 
ergy dipole magnitude in the context of isotropic ACDM 
along with the observed dipole magnitude indicated by 
an arrow. As expected from eq. (|2.15p only 4.75% of 
the simulated isotropic datasets had a dark energy dipole 
magnitude larger than the observed value. This is con- 
sistent with eq. (|2.15|) which indicates that the statisti- 
cal significance of the existence of a dark energy dipole 
is about 2c. In Fig. [5] we show the probability dis- 
tribution of the angular distance of the isotropic simu- 
lated dipoles from the observed fine structure constant 
dipole discussed in section |2A1 Only 6.12% of the Monte 
Carlo datasets had such an angular distance smaller than 
the observed one. The probability for a Monte Carlo 
isotropic Union2 dataset to have both a dipole magni- 
tude larger than the observed one and an angular sepa- 
ration from the fine structure dipole smaller than the ob- 
served one is 0.98%. This is larger than the anticipated 
value of 0.0612 x 0.0475 = 0.29% due to the nonuni- 
form distribution of the Snla in the sky. The conver- 
gence of these probabilities as we increase the number 
of Monte Carlo simulated isotropic datasets is shown in 
Fig. [5] Clearly the number of simulated datasets consid- 
ered (10 4 ) is enough to achieve the convergence of the 
required probabilities. We estimate the combined prob- 
ability that both dipoles have magnitudes larger than 
the observed and angular separation smaller than the ob- 



served in the context of isotropic underlying models to be 
less than 0.01% x 0.98% ~ 0.0001% where the first factor 
comes from the magnitude of the fine structure constant 
dipole estimated in the previous subsection. 

3. REDSHIFT TOMOGRAPHY 

In the previous section we have shown that the dipole 
anisotropy model provides a significantly better fit than 
the isotropic model for both the fine structure constant 
and for the dark energy spatial distributions. We also 
demonstrated that the two dipole directions show a re- 
markable coincidence. In this section we focus on identi- 
fying the redshift ranges in which these effects are more 
prominent. We use two approaches: a redshift bin ap- 
proach and a variable upper redshift cutoff approach. 
In the redshift bin approach, we divide each dataset in 
three redshift bins of approximately equal number of dat- 
apoints and perform an analysis similar to that of the pre- 
vious section in each bin. Thus we compare the results 
of each bin with respect to the quality of data (errorbar 
sizes), the dipole magnitudes and the dipole directions. 
In the variable upper redshift approach we start with 
truncated datasets with an upper redshift cutoff consist- 
ing of about 1/2 of the datapoints. Then we increase the 
upper redshift cutoff in five steps so that in the final step 
the almost full dataset is obtained. We analyse each one 
of the six cumulative dataset parts with respect to their 
dipole magnitudes and their directions. 

In TableUJwe focus on the Keck+VLT sample and show 
the redshift ranges of each redshift bin and of each one 
of the six cumulative redshift parts. For each redshift 
range we show the corresponding best fit monopole mag- 
nitude, the dipole magnitude, the direction of the best fit 
dipole in galactic coordinates and its angular separation 
from dark energy dipole (obtained from the full Union2 
dataset). A similar redshift tomography for the Union2 
dataset is presented in Table [TTI in which we also consider 
a redshift bin that is common to the two datasets (last 
line) . The directions of the best fit dipoles for each one of 
the redshift ranges considered in Tables HI and [TTI is shown 
in Fig. [7] (the cumulative redshift parts are separately 
connected according to increasing redshift cutoff). The 
uncertainties shown in Tables HI and [TTI are obtained using 
the covariance matrix approach. We have checked that 
they are in good agreement with the corresponding la 
errors obtained from the Monte Carlo simulations. 

The following comments can be made based on the 
results shown in Tables Q] and [TTJ 

• The redshift bin with the smallest la errors (best 
data quality) for the Union2 data is the low redshift 
bin {z G [0.015,0.14]). The corresponding best data 
quality redshift bin for the Keck+VLT dataset is 
the high redshift bin (z £ [2,4.1798]). 

• These best quality redshift bins also have the best 
dipole alignment and the most statistically signif- 



6 



(b=90) 




(b=-90) 



FIG. 3: Union2 datapoints and Dark Energy dipole direction. Datapoints in three different redshift bins are represented with 
different shapes. For comparison the direction of the a-dipole obtained from the best fit of the Keck-VLT data is shown with 
a star. The light blue blob represents the 1-a error on the Dark Energy dipole direction. 
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TABLE I: Keck+VLT data: Monopole, dipole magnitude and direction and angular distance from the Dark Energy dipole 
in several redshift ranges. The angular distance with respect to the Union2 dipole is referred to the full redshift case for the 
Union2 dataset. We don't include the common range bin (0.2223 < z < 1.4) since it differs from the fifth line in the table only 
by one datapoint (with z = 0.2223). 



icant deviation of the best fit dipole magnitudes 
from isotropy. It is therefore important to improve 
the quality of data in the other redshift bins in or- 
der to clarify whether the dipole trend is also strong 
in these bins where the data quality is lower. 

The above points are also demonstrated in Figs. [5J 
[9] and [TO] For each one of the redshift bins consid- 
ered, we show in Fig. [5] the fraction of isotropic Union2 
datasets that exceed the observed dipole dark energy 
magnitude and also have a smaller angular distance from 
the Keck+VLT dipole than the actually observed angu- 
lar distance. Clearly, this fraction is significantly lower 
for the lowest redshift bin which implies that the dipole 
behaviour and alignment is most significant for this red- 
shift bin. In Fig. [9] we show the angular separation of 
each Union2 redshift bin from the best fit dipole direction 
of the full Keck+VLT dataset, as a function of redshift 
range for each Union2 bin. Clearly, the lowest redshift 



bin which also has the smallest angular separation er- 
ror is the one that has its dipole best aligned with the 
Keck+VLT dipole. Finally, in Fig.[10]we show the angu- 
lar separation of each Keck+VLT redshift bin from the 
best fit dipole direction of the full Union2 dataset, as a 
function of redshift range for each bin. Clearly, the high- 
est redshift bin has the lowest error and good alignment 
with the Union2 dipole. In this case however, the best 
fit dipole direction appears to be more consistent among 
the three redshift bins while the errorbars for the dipole 
direction are significantly smaller than the Union2 case. 

The above choice of redshift bin ranges has been made 
by demanding approximately equal number of datapoints 
in each one of the three redshift bin. In Ref [15| two 
redshift bins were considered: a high redshift bin with 
z > 1.6 and a low redshift bin with z < 1.6. The moti- 
vation for this redshift division comes from the fact that 
high redshift absorber spectra are dominated by differ- 
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datapoints 


0.015 < z 


< 1-4 


2.0 ± 2.2 


1.3 
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0.6 


-15.1 ± 11.5 


309.4 ± 18.0 


11.3 


± 


17.3 


557 


0.015 < z 


< 0.14 


2.6 ±3.4 


1.7 


± 


0.8 


-10.1 ± 15.1 


308.8 ± 22.8 


11.6 


± 


22.1 


184 


0.14 < 2 


< 0.43 


2.6 ±5.6 


1.2 


± 


1.9 


-10.7 ± 28.7 


291.4 ± 37.2 


28.6 


± 


36.7 


186 


0.43 < z 


< 1.4 


0.7 ±4.3 


0.9 


± 


0.8 


-25.1 ± 30.6 


34.3 ± 75.7 


70.6 


± 


68.7 


187 


0.015 < z 


< 0.23 


3.3 ±2.9 


1.8 


± 


0.7 


-8.5 ± 12.4 


302.2 ± 16.6 


18.3 


± 


16.0 


239 


0.015 < z 


< 0.31 


3.8 ± 2.9 


1.9 


± 


0.7 


-7.6 ± 11.6 


307.0 ± 14.7 


13.9 


± 


13.8 


292 


0.015 < z 


< 0.41 


o.U ± A. 1 


1.8 


± 


0.7 


— 14.4 zh 1U.O 


OHO R -L- 1 A A 

oUo.O ± 14.4 
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14.1 
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-14.9 ± 12.7 


301.3 ± 18.8 


18.9 
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18.2 


406 


0.015 < z 


< 0.64 


2.1 ± 2.4 


1.4 


± 


0.6 


-16.0 ± 11.0 


305.3 ± 16.9 


15.4 


± 


16.2 


464 


0.015 < z 


< 0.89 


2.2 ± 2.3 


1.4 


± 


0.6 


-15.6 ± 10.4 


309.8 ± 16.0 


11.1 


± 


15.3 


519 


0.2223 < 


z < 1.4 


1.2 ± 2.5 


0.51± 0.48 


-44.0 ± 62.5 


59.3 ± 147.6 


88.2 ± 110.6 


319 



TABLE II: Union2 data: Monopole, dipole magnitude and direction and angular distance from the a-dipole in several redshift 
ranges. The angular distance respect to the a-dipole is referred to the full redshift case for the Keck-VLT dataset. 





mK/y(10- 6 ) 


d K/v (W~ 5 ) 


b d K/v (°) 




RA dK/v (hr) 


dec d K/v (°) 


dp 


0.2223 < z < 1.6 
1.6 < z < 4.1798 


-3.9±1.1 
1.1±1.4 


0.57 ± 0.26 
1.39 ± 0.35 


-16.4 ± 15.1 
-10.5 ± 7.8 


336.4 ± 22.9 

325.5 ± 11.9 


18.1 ± 2.0 
16.6 ± 1.4 


-57.3 ± 20.9 
-63.0 ± 10.2 


148 
145 



TABLE III: Monopole, dipole magnitude and direction in the low and high redshift ranges for the Keck±VLT data. The 
results have been obtained fitting the data using the same values of a ran d as in [l5[ (three values) and removing the two outliers 
as identified by [HI]. These results are almost identical with those of Ref. [3] which provides a good test of our analysis. We 
have checked that using a single value a ra nd = 1-0 for all datapoints (as done in the rest of our analysis) leads to consistent 
results and affects mainly the error bars which become somewhat larger. 
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FIG. 4: Distribution of dark Energy dipole magnitudes ob- 
tained from the Monte Carlo simulation. The arrow points to 
the position of the observed best fit value and the light green 
area indicates fraction of the Monte Carlo datasets that give 
a dipole magnitude bigger than the observed best fit one. 
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FIG. 5: Angular distances between the observed a-dipole 
direction and the dipole direction obtained from the Monte 
Carlo simulations on the Union2 data. The arrow points to 
observed angular distance value and the light green area rep- 
resents the Monte Carlo datasets that give an angular distance 
smaller than the observed one. 



ent absorption lines compared to low redshift absorption 
spectra. Thus a test for possible systematics could be to 
divide the whole sample and compare the two resulting 
dipoles. If the two dipoles are consistent with each other 
then this is an indication that no systematic errors are 
hidden in the different absorption lines. No systematic 
errors were found in Ref [15| using this approach. 

In order to test the consistency of our results with 
those of Ref. [l5| we have used the same two bins and 
constructed Table IIIII which is to be compared with the 
corresponding Table 3 of Ref. . In an effort to repro- 
duce the results of [HI we have used the same values of 
fraud and ignored two outlier datapoints. Our results are 
almost identical with those of Ref. [l5| and this provides 
a good test of the validity of our analysis. 



4. PHYSICAL MECHANISM: EXTENDED 
TOPOLOGICAL QUINTESSENCE 

If the observed coincident large dipole anisotropics are 
due to a physical mechanism and not to systematic or sta- 
tistical fluctuations, then it is of particular interest to in- 
vestigate what could be a physical model that would give 
rise simultaneously to these coincident dipoles. Such a 
mechanism could involve for example an inhomogeneous 
scalar field which couples to electromagnetism through a 
non-minimal coupling and whose potential energy could 
provide the dark energy required for accelerating expan- 
sion. Due to negative pressure such a scalar field would 
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FIG. 6: Percentage of Union2 Monte Carlo dataset satisfying 
different constraints as a function of the number of Monte 
Carlo datasets considered. The points labelled as "magni- 
tude" represent the fraction of Monte Carlo datasets that 
give a dipole magnitude larger than the observed one. Those 
labelled as "ang. distance" represent the fraction of Union2 
Monte Carlo datasets that have an angular distance from 
the observed a-dipole smaller than the observed angular dis- 
tance. The label "both" refers to the fraction of Monte Carlo 
datasets that satisfy both the previous constraints. With 
"product" we label the points that represent the product of 
the first two percentages. 



tend to quickly become homogeneous and isotropic on 
Hubble scales. However, nontrivial topology would nat- 
urally generate sustainable inhomogeneity [25| of such a 
scalar field. 

For a proper potential, the scale of the inhomogene- 
ity would be the observationally required Hubble scale. 
In such a Hubble scale topological defect, an off cen- 
tre observer would observe aligned dipoles in both dark 
energy and the fine structure constant. For a large 
enough core scale, such a defect would become effectively 
homogeneous and indistinguishable from ACDM . The 
dipole nature of observations of off-centre observers lo- 
cated in spherically symmetric inhomogeneities has been 
discussed in detail in Refs [H, [27} • 

In the case of no coupling to electromagnetism, this 
mechanism was studied in detail in Ref. 25] (topolog- 
ical quintessence). Topological quintessence is an ex- 
tension of the well known corresponding inflationary 
model: topological inflation 28j. In what follows we 
present some qualitative features of the extended topolog- 
ical quintessence model and we postpone a more detailed 
study for a later publication. 

Consider the action 



1 M*R- \{d^ a Y-V{®)- 
1 



(4.1) 



where M p 2 = 8irG is the reduced Planck mass, C m is 
the Lagrangian density of matter fields, $ a (a = 1,2,3) 



is an 0(3) symmetric scalar field, B(^>) is a non-minimal 
coupling to electromagnetism and 



1 



V{$) = -A($ 2 - ?y 2 ) 2 , $ = V$ Q $ a 



(4.2) 



We assume the existence of a Hubble scale global 
monopole formed during a recent cosmological phase 
transition. The vacuum energy density in the monopole 
core and the size of the core are determined by the 
two parameters of the model rj (the vacuum expecta- 
tion value) and A (the coupling constant). The global 
monopole field configuration is described by the hedge- 
hog ansatz 



$ a = $( r ,t)r a = $(r, t) (sin 9 cost/?, sin 9 sin tp, cos ( 



shown in Fig. 11(a) with boundary conditions 
$(0,i)=0 , $(oo,i) = ?7, 



(4.3) 



(4.4) 



where rj is the scale of symmetry breaking assumed to be 
such that (H 



Xrf_ 
3H 2 



1 



(4.5) 



In eq. (|4.3|) we have allowed for a time dependence hav- 
ing in mind a cosmological setup of an expanding back- 
ground. For a slowly evolving global monopole configu- 
ration, the size of the core is approximated by 

^-A-^V 1 , (4.6) 
while the vacuum energy density in this core region is 

P c ° re * ■ (4.7) 

For a core size much larger than the Hubble scale, the 
model reduces to ACDM . Therefore the constraints im- 
posed on inhomogeneous matter models (29j are not appli- 
cable to this class of inhomogeneous dark energy models. 

The general spherically symmetric spacetime around 
a global monopole may be described by a metric of the 
form 

ds 2 = -dt 2 + A 2 {r,t)dr 2 + B 2 {r,t)r 2 {d9 2 + sin 2 9d<p 2 ). 

(4.8) 

A detailed analysis of the cosmological evolution of the 
above metric in the presence of the global monopole and 
matter is presented in Ref. [25| with £?($) = 1 (see also 
Ref. [HQ). 

Consider now a non-minimal coupling of the form 



B($) 



1 



$2 



(4.9) 



where £ is constant. The fine structure 'constant' is re- 
lated to the coupling B(&) as 



a($) 



47rB($) 2 



(4.10) 
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FIG. 7: Fine structure a and Dark Energy dipole directions for the different redshift bins. The 'stars' denote the bins 
corresponding to the redshift range that is common to the Keck+VLT and Union2 samples. For this range however, the dipole 
uncertainty obtained from the Union2 data is very large (see last line of Table (tTTjl . 
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FIG. 8: Percentage of the Monte Carlo Union2 datasets that 
give both a Dark Energy dipole larger than the observed value 
and an angular distance between the Dark Energy dipole and 
the Q-dipole smaller than the observed angular distance. The 
result is plotted for the full and partial redshift bins. 
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FIG. 9: Angular distances (with errors) between the Dark En- 
ergy dipole obtained in the different redshift bins and the a- 
dipole obtained from the full redshift range Keck+VLT data. 
Notice that the lowest Union2 redshift bin dipole has the best 
alignment with the Keck+VLT dipole and also has the small- 
est error. 



where eo is the bare charge that remains constant 
throughout the cosmological evolution. Therefore for 
small values of $/ry we have 



a J " rf 



(4.11) 



where $o is the field magnitude at the location of the ob- 
server. The dipole directions shown in Figs.[T]and[3]corrc- 
spond to higher value of a and lower accelerating expan- 
sion (brighter Snla compared to ACDM ) respectively. 
Thus, in the extended topological quintessence picture, 
for an off-centre observer, this would be the direction 



pointing away from the global monopole core where the 
potential energy of the monopole is lower and the field 
magnitude $ is larger. In order to have a higher value of 
the fine structure constant in the same direction we need 
£>0. 

In Fig. [IT] we illustrate the location of an off-centre ob- 



server with respect to the monopole core. In Fig. 11(a) 



we plot the observer location along with the field magni- 
tude and direction denoted by the arrows at each point 
of the x — y plane. Clearly, the field magnitude is smaller 
towards the centre of the monopole and this justifies the 
variation of a in that direction. Similarly, in Fig. 11(b) 



we show the energy density distribution of the global 
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FIG. 10: Angular distances (with errors) between the a-dipole 
obtained in the different redshift bins and the Dark Energy 
dipole obtained from the full redshift range Union2 data. No- 
tice that the alignment of all Keck+VLT redshift bin dipoles 
with the full Union2 dipole is consistent with each other and 
similar to the alignment of the full Keck+VLT dataset. 



monopole and the location of the observer. Clearly, there 
is higher dark energy density towards the monopole cen- 
tre and this justifies the higher acceleration rate in that 
direction. 

We postpone a detailed reconstruction of the global 
monopole potential V(&) and coupling -B($) for a later 
publication. A comparison of the quality of fit for dif- 
ferent topological defect geometries could also be made 
(global vortex or thick domain wall). We also stress that 
our extended topological quintessence approach is dis- 
tinct from the model of Refs. [13, Hl[ where thin domain 
walls were considered in an effort to explain the spatial 
variation of a. The approach of Refs. [2(J Hl[ does not 
address the dark energy dipole and predicts a non-dipolc 
anisotropy for a. The dipole anisotropy however has been 
shown [30( to provide a better fit to the Keck+VLT data 

The double wall model 



that the single wall model [21 1 

[2fj| however, involving three additional parameters, has 
been shown to provide a better fit than the dipole model. 

The extended topological quintessence monopole dis- 
cussed above is also distinct from the varying a defects 
p^H based on Bckensteins's theory [32]. According to 
this model, the electric charge e (and therefore also the 
fine structure constant a) is promoted to a dimension- 
less scalar field tp ~ In e with zero potential and a ki- 
netic term multiplied by a large dimensionful parameter 
oj. This is similar to the corresponding extension of Gen- 
eral Relativity along the lines of the Brans-Dicke theory 
where Newton's constant is promoted to a scalar field. 
As in the Brans-Dicke theory, the parameter uj is used to 



partly freeze the dynamics of <p so that the charge varia- 
tion in spacetime becomes consistent with observational 
and experimental constraints [H, HH . The dynamics of 
the charge field ip affects the dynamics of the gauge field 
which in turn affects the dynamics of any scalar field 
<& that couples to via a gauge symmetry. Varying 
a defects [3l[ are gauged defects formed if the vacuum 
manifold of <f> has a non-trivial homotopy group and their 
dynamics is indirectly affected by the dynamics of the 
fine structure constant (and of tp). A potential source 
of severe constraints for this class of defects is that they 
predict massive photons (spontaneous breaking of elec- 
tromagnetism) in regions away from the defect core. 

In contrast to these varying a defects, in extended 
topological quintessence, the defect is global and is 
formed by the same field that represents a. These 
are global defects non-minimally coupled to electromag- 
netism. As a simple example in Minkowski spacetime, 
consider a global vortex non-minimally coupled to elec- 
tromagnetism. The dynamics of the complex scalar field 

is determined by the Lagrangian density 

C = (0„ $)* (8" *) - j Ffj, u - V($) , (4.12) 

The field equations obtained by variation of <J>* and 
are 

W = -^-i«F M „i- (4.13) 



and 



0„ =0, (4.14) 

For a non-minimally coupled global vortex to form we set 



and 



|$|2 



(4.15) 



(4.16) 



We now use the global vortex the ansatz allowing for a 
coaxial magnetic field 



$ = /(r)e' 
A g = a(r) , 



(4.17) 
(4.18) 



where f(r) and a(r) are real functions of r and all other 
components of a M are set to zero. We thus obtain the 
static field equations for f(r) and a{r) as 
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FIG. 11: (a) The observer location (thick dot) along with the field magnitude and direction denoted by the arrows at each point 
of the x — y plane. Clearly, the field magnitude is smaller towards the centre of the monopole and this justifies the variation of 
a in that direction, (b) The energy density (p) distribution of the global monopole and the location of the observer. Clearly, 
there is higher dark energy density towards the monopole centre and this justifies the higher acceleration rate in that direction. 



1jL( ^€\-( r f-- r ^ x i A 

r dr \ dr J 
dr \ r dr 



ldB{f 2 ) (1 d 



df 



dr 



(ra) 



(4.19) 
(4.20) 



since 



F» u = 2F re F r , 



1 d 

r dr 



(ra) 



(4.21) 



The corresponding energy density of the vortex is: 



P 



1 

^2 



d (r a) 
dr 



(4.22) 

If there is no external source of electromagnetic field and 
if B(f 2 ) > everywhere, we obtain the usual global vor- 
tex solution f(r) = fo(r), ol(t) = 0. However, if there are 
regions of space where B(f 2 ) < 0, an instability develops 
which proceeds with spontaneous creation of electromag- 
netic field in the region where B(f 2 ) < 0. For example 
for B(f 2 ) = 1 — qV(f 2 )/ri 4 , an instability develops in the 
core, for large enough values of q. 

If there is an external source of electromagnetic fields 
(e.g. a localised magnetic field in the z direction), then 
the profile of f(r) will be affected in accordance with 



eq. (|4.19l) and a local additional variation of a will oc- 
cur. Thus, a robust prediction of this class of models is 
a correlation between regions of strong electromagnetic 
fields and variation of a. The non-observation of such 
variation could impose strong constraints on the form of 
the coupling B(f 2 ). The detailed investigation of these 
constraints and their consistency with the form of B(f 2 ) 
required to explain the observed a dipole represents an 
interesting extension of this project. 



5. CONCLUSIONS 

We have used the Keck+VLT dataset and the Union2 
dataset to show that the value of the fine structure con- 
stant and the rate of accelerating expansion are better 
described by coinciding dipoles than by isotropic cos- 
mological models. The key feature of our analysis is 
that it applies identical method (fit to dipole+monopole 
anisotropy) to both the Keck+VLT dataset and the 
Union2 dataset. This consistency, combined with the 
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apparent dipole nature of the anisotropy, has allowed a 
consistent comparison of the two dipoles. 

Using Monte Carlo simulations and covariance matrix 
error estimates, we find that the probability that these 
coinciding dipoles are both produced in the context of 
a cosmological model where fine structure constant and 
dark energy are isotropic and uncorrelated is less than 
one part in 10 6 . A redshift tomography analysis di- 
viding the two datasets in three redshift bins revealed 
that the highest data quality redshift bins correspond to 
low redshifts for the Union2 sample and high redshift 
for the Keck+VLT data. The dipole direction for the 
Keck+VLT data depends weakly on redshift while the 
Union2 dipole direction depends more strongly on red- 
shift and it is the low redshift (and lowest error) bin that 
is best aligned with the Keck+VLT dipole. The direc- 
tional uncertainty is significantly larger for the medium 
and higher redshift Union2 dipoles. It is therefore im- 
portant to improve the quality of intermediate and high 
redshift Snla data in order to further test the alignment 
of the dark energy dipole with the fine structure constant 
dipole. 

An important issue that we have not addressed in the 
present paper is the effect of systematic errors of the 
Keck+VLT sample. This issue has been addressed in 
detail in Ref. [15[ where no significant source of system- 
atic errors was identified. The main concern has been 
the possibility of careless merging of the two datasets 
(VLT and Keck) which in principle have different sys- 
tematics and effectively cover opposite hemispheres of 
the sky which coincide with the direction of the identi- 
fied Keck+VLT dipole. The concern therefore is that the 
large identified dipole magnitude originates from a hid- 
den difference in systematic errors between the VLT and 
Keck samples J22J. According to Ref. this does not 
appear to be the case for the following reasons: 

• the dipole directions at high and low redshifts are 
in agreement (this is confirmed in our study too as 
shown in Figs. 171 and [T0|) ; 

• the directions of the dipoles fitted by the VLT and 
by the Keck samples separately are in agreement; 

• the absorbers that are common to both the VLT 
and the Keck sample provide consistent values for 



logical quintessence). In such a model, an off-centre ob- 
server with respect to the Hubble scale core of a global 
monopole would naturally observe large aligned dipoles 
for the fine structure constant and dark energy. In fact 
it should be possible to reconstruct both the scalar field 
potential and the non-minimal coupling form using the 
Keck+VLT and the Union2 samples. 

A robust prediction of the non-minimally coupled de- 
fect model is the weak dependence of the value of a on 
the existence of local strong magnetic fields as discussed 
at the end of the previous section. Another interesting 
prediction is the existence of peculiar velocities in the di- 
rection away from the center of the global monopole due 
to the repulsive effects of antigravity (negative pressure) 
in the defect center. An off-center observer would expe- 
rience this Hubble scale flow as a dipole dark flow. Such 
dipole dark flow has indeed been observed 0, H[ and it 
is attributed to the existence of a Great Attractor which 
could be present on Gpc scales (perhaps even at a neigh- 
boring universe 35] ) . In our model such a dark flow could 
be due to a 'Great Repulser' whose role would be played 
by the core of the Hubble scale non-minimally coupled 
defect. The predicted direction of such a flow should be 
away from the defect core (Great Repulser) in the di- 
rection of maximum deceleration {b = —15.1° ± 11.5°, 
I = 309.4° ± 18.0°) (see Table HJ). The direction of the 
observed dark flow is (6 = 8° ±6°, I = 287° ±9°) 2] which 
is consistent within la with the direction of the dark en- 
ergy and a dipoles. It also points towards the region of 
lower acceleration as predicted by our model. A robust 
prediction of our model with respect to the dark flow is 
that it should reverse direction at large enough redshifts 
as we start seeing on the 'other side' of the 'Great Re- 
pulser' (defect core). 

A detailed investigation of the consistency of the above 
predictions with cosmological observations is an interest- 
ing extension of the present analysis. 

Numerical Analysis Files: The data, Math- 
ematica and C-l — hprogram files used for the nu- 
merical analysis files may be downloaded from 
http : / /leandros .physics . uoi . gr/def sdipoles. 
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